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Abstract. Atomic response to a probe beam can be tailored, by creating coherences between 
atomic levels with help of another beam. Changing parameters of the control beam will change 
the nature of coherences and hence the nature of atomic response as well. Such change can 
'jrt ' depend upon intensity of both probe and control beams, in a nonlinear fashion. We present 

^ [ a situation where this nonlinearity in dependence can be precisely controlled, as to obtain 

O ' different variations as desired. We also present a detailed analysis of how this nonlinear 

c/3 . dependency arises and show that this is an interesting effect of several Coherent Population 

^^1 Trap(CPT) states that exist and a competition among them to trap atomic population in. 

^: 
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^ ! Response of a material medium to an electromagnetic field usually exhibits a linear 

^ ! dependence on the intensity, at low intensities. When the intensity is sufficiently high, 

■*>-> ! intensity dependant components of refractive index and absorption become dominant and the 

O ! response of the material starts becoming increasingly nonlinear. While this is usually true for 

^! bulk medium, dilute atomic gases show a nonlinearity even at lower intensities. This is due to 

Jh '. the fact that a near resonant laser prepares the atom in a coherent superposition of its energy 

levels. The response of an atom in such a superposition is not merely a sum of responses 
from the atom in its constituent states and can be dramatically different. Thus, by creating a 
superposition state, a control laser can alter the behaviour of the atom towards a probe laser 
5^ \ propagating through such a medium. Such coherence induced control of an atomic response 

can also be interpreted as a coupling between photons of the control beam and those of the 
probe beam, giving us the same picture as in traditional nonlinear optics |II]|,|l2|. Although 
a nonlinear response at higher intensities is also observed, due to anharmonic interaction 
of the atom with the laser, the coherence induced nonlinearity is much more efficient since 
it happens at much lower intensities. This phenomena has found several applications such 
as achieving giant nonlinearities ||3J, efficient generation of of Vacuum Ultra Violet (VUV) 
using third harmonic[4|, suppression of nonlinear responses [5|, light induced focusing |l6| 
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and many more. A fast variation of atomic response between a complete transparency to a 
complete opacity can also be achieved which would effectively make an optical switch using 
the medium f 71 . 

A recent study [|8j| experimentally investigates several A-like systems for dependency 
of atomic response on probe intensity and shows existence of a threshold intensity where 
linearity breaks down. On the other hand, we present a slightly different configuration in 
this paper which exhibits an almost completely nonlinear response. Further, this nonlinearity 
is precisely controllable, by manipulating the intensity and frequency of the pump beam. 
Our configuration consists of two hyperfine states with quantum number F=l and its Zeeman 
manifolds, totally amounting to six levels. Such a system is often encountered in atoms such 
as ^^Rb or Sodium, which is used in many atom-laser experiments. In other words, this 
is a familiar system which exhibits a sufficiently rich nonlinear response under appropriate 
conditions, as explained in this paper. We also present in detail the algorithm used to obtain 
numerical solutions which describe the dynamics of this system. Due to involvement of 
two lasers of slightly differing frequencies, the most used method of transformation into a 
rotational frame will not eliminate explicit time-dependencies in the interaction Hamiltonian 
and hence we resort to a continued fraction approach to solve them. This method, being non- 
perturbative in approach is valid for any intensities of pump and probe beams. Thus, we are 
not confined to regimes of strong pump and weak probe alone, and are able to continuously 
investigate the behaviour of the atom as probe beam's intensity is increased from almost zero 
value to a value larger than that of the pump beam. Augmenting the study with an analysis 
using Semi-classical Dressed states, we are able to show that the complex response shown by 
the atom is due to existence of many Coherent Population Trap (CPT) states which compete 
with each other for atomic occupation. 

CPT states are formed when superposition of atoms are prepared in such a superposition 
of its bare states, that the probability of transitions from each of its constituent states to a 
common external state cancel each other. Thus, the atom can not absorb a photon to make a 
transition to the common external state is thus trapped |J9|- However, there may exist other 
states to which transitions not canceled, and the CPT state can be depopulated by using the 
appropriate laser. Our configuration, due to existence of many degenerate levels, leads to 
formation of more than one CPT state. As the parameters of the lasers, such as frequency 
and intensity are scanned, the atoms are transferred from one CPT state to other. The rates 
of optically pumping the population to a particular CPT state or out of it depends upon the 
nature of the state. Therefore, in the intermediate states, a competition between different CPT 
states manifests, as to trapping the population. This phenomena is observed as a variation in 
probe absorption and dispersion, which varies nonlinearly with variation of intensities of the 
beam, as shown in this paper. Furthermore, the nature of nonlinear response shown to one of 
the beams can be modified by varying the detuning or intensity of the other beam involved. 
Thus, a precise control over the nonlinearity is possible, which may have several applications 
in information processing using nonlinear properties of the medium. 

This communication is divided into three parts. We describe the atomic configuration 
in detail in the following section, and also describe the method we have adopted to solve 
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the density matrix equations. In the subsequent section, we present the results obtained under 
different conditions of pump and probe intensities. From this numerical results, we attempt an 
insight into the phenomena that occurs. We obtain the semiclassical dressed states in section 
IV and substantiate our conclusion of competition between different CPT states. 

1. The configuration 

The atomic level configuration consists of six levels as shown in figure 1 . The levels labeled 
\sq): \s+): \S-) form the triplet of the ground state Fg=l state with magnetic quantum 
numbers mp = 0, ± 1 respectively. Those, that are labeled |eo), k+) and |e_) are similarly the 
Zeeman sublevels of the excited state F^ = 1 , with their magnetic quantum numbers mp = 0,1 
and —1. For sake of convenience the levels Fg=l and Fg = 1 could be hyperfine levels of the 
ground state 5S W2 and of the first excited state SPj /2, respectively of ^''Rb. Any two hyperfine 
states with F=l, of any other alkali atom, such as Sodium for instance, will give identical 
results. 

Two beams, plane polarized in orthogonal directions are incident on this system. For the 
probe beam, we use the one which is plane polarized along the z-axis of quantisation. This 
beam couples the transitions |g_) ^^ |e_) and |g_|_) ^^ |e_|_). The dipole transition between 
I go) ^^ |eo) is forbidden due to vanishing Clebsch-Gordon coefficients. The other beam is 
plane polarized in xy plane, thus forming a combination of two circularly polarized lights with 
helicities a+ and o . This beam therefore couples transitions |g±) ^^ \qq) and |go) ^-^ |e±). 
We refer to this as the pump beam. We then analyse the response of the medium to the probe 
beam, as modified by the pump beam. 

In the semi-classical approach, the pump and the probe beam are represented by their 
electromagnetic fields Ei = •Eiexp{i(Dit — ki.z) + c.c. and Ep = 'Epexp{i(x)pt — kp.y) + c.c. 
respectively, where £/ is the slowly varying envelop of the pump field and 'Ep, that of the 
probe field. (0/ and odp are their respective angular frequencies. The pump beam is assumed 
to propagate in z direction, while probe beam can propagate either in x or y directions to give 
the required polarization. We assume it to be y direction for the sake of convenience. For rest 
of this paper, we refer these two beams by their respective Rabi frequencies Q/ p = dtj.'Eip 
where i = \eo,±) and j = |^o,±)- 

A closely resembling configuration has been analyzed and reported by Petrosyan et. al., 
Q showing an efficient two-photon process occurs in the media, with one photon from the 
pump beam and another from the probe. The authors present study of propagation effects 
in a situation where the \g±) and \e±) levels are Zeeman shifted. Our analysis differs from 
theirs in three aspects. Firstly we look at the steady state phenomena, which can model 
most experiments involving cw lasers as well as very long pulses. Secondly we analyze a 
nearly-degenerate system with no Zeeman shifts involved, which, as we show, exhibits a wide 
variety of non-linear phenomena. Thirdly, we derive a general method for obtaining numerical 
solutions of such systems with no restrictions on intensities of pump and probe beams. The 
method obtains results to all orders of magnitude in pump and probe beams and hence can be 
adopted even when probe and pump beams are of comparable intensities, or even when probe 
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is stronger than the pump. We present the procedure involved to obtain numerical solutions 
below. 

The Hamiltonian governing the system is given by 

^=^Q + 9{int, (1) 

where, 

:?/o = a)e|e+) (e+l + (Oe|e_) (e_ I + (Oeko) (eol 

+ ^s\g+){g+\+^g\g^){g-\+^g\go){gQl (2) 



-^pe-^^"^+^'^py{\e+){g+\-\e^){gA)+H.C.. (3) 

COg^e are the eigen-energies of the states |go,±) and |eo,±) respectively in terms of angular 
frequencies and h=\. The coupling factors are defined as 

fora= |eo,±) and Z? = |go,±) 

Relative sign difference for transition between |e_|_) ^^ 1^+) and |e-) ^^ |^-) in Q arises 
due to a sign difference between dipole matrix elements such that 

ag g agj^gj_. 

The evolution of 6x6 density matrix of the system is governed by the Liouville equation 

dtp = -i[!}{,p]+Lp. (4) 

ip denote the terms involving decay terms that originate from both spontaneous 
emission and coUisional mixing. We assume that collisions do not bring about any significant 
mixing between \e) and |^) states although coUisional transfer within the Zeeman substates 
of \e) and \g) are significant. That is, most collisions are of the spin-flip kind. The equations 
for each term of density matrix p are explicitly written in Appendix. Traditional method 
of obtaining solutions for the equations involve a transformation into a rotating frame of 
reference such that the right hand side of equation Q has constant coefficients. However 
using either of the two frequencies involved, namely that of probe or the pump beam, have 
oscillating coefficients. In other words, there exist two rotating frames of reference with 
rotational frequencies of (0/ and cOp, and choosing either one of them as stationary will result 
in the other one rotate at a frequency equal to the difference between the two of them. 

To solve such a system, we chose a rotational transformation with the pump frequency co; 
such that the transitions connected by pump beam is seen stationary §. These are the transitions 

§ Choosing the system which is resonant with co^ would equally serve. 
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|^_) ^^ |eo) ^^ 1.?+) and |e_) ^-> |^o) ^^ k+)5 which are the three levels configurations A and 
V, which are independently studied in the past. The probe beam couples these two subsystems 
whose transitions are now investigated. 

Accordingly, a rotational transformation is operated such as 

pij -^ pijcxp{~i(i)it); for / 7^ j 

Pa -^ Pa 

for ij= 1,...,6 

The equations in such a rotated frame become 

dtp=Ap + e'^""-''p^'M+p + e-''^'^'-"'"^'M^p + Io. (5) 

A is a 35 X 35 matrix that involves all interaction terms which are due to (0/ field alone, 
and those related to free evolution. These terms are seen as stationary in the rotating frame 
chosen. On the other hand, terms in M_|_ and M_, which are also 35 x 35 matrices contain 
terms coupled to Q.p field and thus have terms with (0/ — cOp in the exponent. Terms with a 
positive exponent (i.e., rotation with a frequency of +((0/ — (Op) are grouped in M_|_ and those 
with a rotation in the opposite sense are grouped in M_ . These terms are go to zero when Q.p 
is zero. 

Iq is the single column matrix of dimension (35x1). The column Iq depends on the way 
the constraint 

Pee + Pgg=l- (6) 

is used to write ©. In © the element pgygo has been eliminated using ®. Explicit 
values for A, M-^, M and Iq that are used in our calculations are given in appendix. 
The solution of © is obtained by assuming 

-|-oo 

p= £ g'"(a),-co„yp(n)_ (7) 

Substituting ©in © and equating synchronous terms on both sides of ©, one obtains 

a^pW = [A-/n(a)/-03p)I]pW+M+p("-i)+M_p("+i)+/o5„o. (8) 

I is the 35 X 35 identity matrix. 

These equations determine the dynamics of the different subsystems that exist in our 
main system. The transfer of populations and coherences from one subsystem to the other can 
be determined from here. In this paper we present results in steady state. The steady state is 
described as the state that is created as a result of equilibrium reached due to the continuous 
exchange of populations and coherences among different subsystems due to the incoherent 
processes like collisions and spontaneous emission as well as coherent interactions due to the 
two fields. Such a steady state is characterized mathematically by setting on left hand side 
^tp^"' = 0- Using this condition leads to three term algebraic recursion relations [101 . which 
can be solved in the following way: 
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We define square matrices X„, X_„, sucli tliat 

pW=X„p("-i); n>0 (9) 

p(-™)=X_„p(-'"+i); m>0 (10) 

and use tliese in the recursion relations for different n, m. One finds 

^n = -r ^ , (11) 

A-in{&i-&p)l + M^Xn+i 

^-m = 7 ^ , (12) 

A + i m{(i)i - (Op)I + M+X_„,_i 
then p^'^-' is given by 

p(0) = ^ . (13) 

This completes consistent mathematical procedure used for the numerical study for 
variety of cases in the six level model. The numerical procedure involves truncating the 
recursion relation dTTT) and (IT2b for sufficiently large values of integers n and m such that the 
result of (IT3t converge to stable values. These are then used to determine other observables of 
the system using ^ and (flOb above. This method, which is non-perturbative in approach, 
can be used for arbitrary intensities of pump and probe beams, as long as result of (fT3t 
converges. The methodology adopted to solve this system can be extended to study several 
other configurations which involve two or more lasers and leading to time dependent terms 
such as those found in this case. We could use the same mathematics to solve a more 
complicated systems where the pump laser interacted resonantly with some transitions and 
off-resonantly with others, which were also coupled by the probe beam. We got consistent 
solutions for that system as well. 

2. Results of numerical solutions 

Using the method discussed in previous section, we numerically solve equation (fT3t and 
obtain total susceptibility as 

x(^)K,£.) = ^^i^-^^- (14) 

Since absorption and refractive index are related to yy^ as inverse Beer length a = 
27r((o/c)Imx(^) and n = 1+271 Rex*-^'' respectively, we compute real and imaginary parts 
of %(^) to qualitatively study the atomic response. We present real and imaginary parts of x^^) 
as a function of probe Rabi frequency ^p, for different values of pump Rabi frequency ^/. 
We notice a pronounced nonlinear variation of both absorption and dispersion of the probe 
beam. For the graphs shown in figure 2, the pump beam is assumed to be resonant with 
the 1^) ^^ \e) transition, while the probe beam is detuned by an amount which is thrice the 
transition linewidth y^g. We do not yet understand why this precise value shows a maximum 
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nonlinearity. Perhaps the cause lies in the nature of different coherent superpositions that are 
formed for different values of co^ and CO/. 

The curves in figure 2, which are labeled 1 to 5, correspond to different values of 
a, = 10"^ 1.0, 2.0,3.0 and 4.0. Parts (a) and (c) show Imx^^), while parts (b) and (d) 
show corresponding variation in Re%(^^ In all of them, Q.p values on the ordinate range from 
a 10^^ to 5.0, i.e., we investigate situations as Qp ranges from a near-zero value, increasing 
upto becoming equal to Q/ and then increasing further. All values used for computation are 
normalized to the coherence decay ygg. Note that this susceptibility is calculated to first order 
in probe beam 'Ep but to all orders in pump beam amplitude £/. 

Note that each curve begins with a different value of %'-^-', indicating that when Q.p is 
nearly zero, the value of x*-^' depends only on Q.i. 

As Q.p is increased, value of %''^-' is determined by instantaneous values of both Q/ and 
Q.p. For very large values of of Q.p, all curves asymptotically reach a single value of X" \ 
irrespective of the Q/ value. In other words, when the two Rabi frequencies differ by a large 
amount, the larger of the two determine the behaviour oix^^> . Varying one of them with the 
other one fixed will show a behaviour of X"^' as seen in curves of figure 2, indicating the 
competition between two Rabi frequencies which control the behaviour. 

Curve number 1 of figure 2a corresponds to f2/ = 10^^, which is effectively a situation 
where only the probe laser is present. Appropriately, this shows a saturation of probe 
absorption as Q.p is increased. Curves 2 and 3, with ^/ = 1.0 and 2.0 respectively, 
show an increased absorption, which eventually saturates. Surprisingly, increasing Q/ still 
further leads to reduced absorption, as shown by curves 4 and 5 of the same figure. The 
accompanying dispersion curves from figure 2 b show a corresponding behaviour. All the 
5 curves asymptotically reach a saturation value of Re X"^' ~ 0. But curve 1 shows an 
early saturation behaviour as compared to others. Curves 2 and 3 begin with an increased 
dispersion and eventually reach the saturation value, although at a much higher value of Q.p. 
Increasing Q/ further, as shown by Curves 4 and 5, dispersion flips the sign and becomes 
negative throughout, eventually approaching zero asymptotically. This means that dispersion 
can be varied between negative and positive values by controlling Q/, and thus a propagating 
probe pulse can be compressed or expanded dynamically. Note that the absolute magnitude 
of Rex*-^-* remains less than one, even when negative. Therefore the refractive index remains 
positive for all values of Q.p. However, we can control the refractive index between values 
less than one and greater than one, by choosing appropriate value of Q-i. 

Plots of figure 2c and 2d show the data for same conditions except that a coUisional 
decay rate of Yc = 0.1 has been introduced. It is clearly seen that X reaches its saturation 
value for a higher value of Q.p, than in case of no collisions. The numerical value of X^^^^ also 
significantly reduces, than its counterpart of no collision case. The qualitative behaviour of 
the variation, including the changes in sign of Re%(^), remains same. It has been shown that 
coUisional processes destroy coherences formed and the laser that creates the coherence is 
required to be stronger in order to balance this loss [fTTl. The effect seen here is similar, since 
it requires a stronger Q.p to saturate x*-^-* when collisions are present. This is an indication that 
the phenomena is indeed driven by atomic coherence, rather than mere optical pumping. 
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We now look at the 'Dressed states' to understand the behaviour seen in above graphs. 
We can safely assume that whenever one of the two lasers are stronger than the other, it is 
the stronger laser which 'dresses' the energy levels and the weaker laser probes the dressed 
atom. Therefore, when Qp << Q/, we look at the atom as it is dressed by fl/ laser. When 
Q.p » Q.I, we look at the dressed states created by Q.p as probed by Q/ laser. Complications 
arise only when both lasers are of comparable intensity when dressed states formed by both 
lasers together has to be taken into account. Seen in this viewpoint, the graphs in figure 2 
indicate the behaviour of the system as it traverses between different set of dressed states. In 
the beginning, when Q.p is small, the system is in dressed states of ^/. As Q.p is increased, the 
system goes through a complex set of dressed states formed due to both Q.p and ^/ lasers and 
finally, reaches the dressed states of Q.p laser. Intuitively, one can understand that for small 
values of Q.p, different values of Q.i form dressed eigenstates with different energies and hence 
the probe absorption is different. When Q.p reaches sufficiently high value, the dressed states 
formed by Q.p dominates and all atoms reach the CPT state formed by Q.p. Thus the atoms 
are trapped and hence both absorption and dispersion reach a value of zero. 

In the following section, we look at the nature of the dressed states to substantiate this 
conclusion. 

3. Analysis in Dressed atom picture 

We consider three separate cases, viz, (i)f2/ >> Q.p, (ii)^/ ~ Q.p and (iii)f2/ << Q.p, and 
analyze the dressed states in each of them. 

(i): ^/ >> Q.p: Dressed states due to Q/: We compute semiclassical dressed states following 
lITlt . The complete Hamiltonian ^ is diagonalized and the diagonalizing matrix is multiplied 
with the Bare state vector to obtain the dressed states as a superposition of bare states. Further, 
this configuration can be divided into two subconfigurations, viz., a A system formed by 
levels \g-) ^^ \eo) ^^ |^+) and a V system formed by |e_) ^^ |^o) ^^ k+)- Subsequently the 
eigensystems for both of them are given as 
forV(e+,^o,e_) 



i: =A 



^o 



II 



A/ipv 



2 

Vo =^{\e+)-\e-)) 

ia+\e+) + a^\e^)±P±\go) 



v^ 



where 



a^+Pl 



a^ =\a+\^ + \a- 



and P± = ^^^. (15) 
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Similarly for A(g+, eo,g-) 
K =0 



1 A/±|3a 



2 
¥^ =^{\8+)-\S-)) 

, _ a*_\s+ )-^l\s-)±Q±\eo) 

Vi — 






a^ = |^_|2 + |^ |2 



Pa =v/A2 + 4a2 



^A 

and e± = ^. (16) 

Also, ^-1- = Q. = ^/. The levels \|/q ^ and \|/q ^ form a ladder of eigenstates, separated 
by a value equal to (O/. There are nine possible coupling between these levels brought about 
by the probe laser, but only four of them have nonzero transition probabilities. Quantum 
interferences causes cancellation of transition in other five cases, as explained by (llTt below 

(\|r^ I <i. Ep I \|/)^ ); cancellation, 

{\\f^\d.'Ep\\\f^); summation, 

(vj;^ I J. E^*^ |\|/^) ; summation, 

(v/± I J. Ep I \/^); cancellation. (17) 

The level scheme, along with allowed transition is shown in figure 3. It is evident from 
( I15I16I ). that the eigenvalues Xq^ and ?Iq ^, and hence the energy differences between them, 
depend upon value Q.± = Q.i. In other words, different Q.i values will induce different energy 
shifts to the relevant energy levels and as a consequence, a probe beam coupling them will 
see different detunings. This will produce different values of %(^) as seen in figure 2, in the 
regions of small Q.p. 

(ii): Q.I « Q-pi Dressed states due to Q.p: The situation of Q.p being stronger than Q.i is 
similar to the one examined by Petrosyan and coworkers QH. In absence of ^/ field, the 
dressed state eigensystem is related to two two-level systems viz., \g±) ^^ |e±), given by 



++) = 


e+) + 


\8+) 


+-) = 


e+)- 


\g+) 


■-+) = 


e-) + 


\g-) 


■— ) = 


e-)- 


\g-) 



(18) 

The Q.p field does not couple |go) ^^ ko) transition. Atoms that reach |go) due to 
spontaneous emission from \e±) are therefore uncoupled from any radiation and are thus 

II The authors of this paper investigate transient coupling between photons of ii/ and £lp through the atom and 
show that this coupling strongly depends upon the Zeeman shift of |^±) and \e±) states by an external magnetic 
field. We present a steady state, zero-field situation which is one of the limiting cases of our investigation 
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trapped in \gQ). However, such atoms are pumped out by Q/ laser, which couples |go) ^^ k±) 
transitions. These transitions are shown by solid lines of figure 4. As Q.p is increased, the 
process of optically pumping the atom into the trap state |^o) dominates over the process of 
pumping it out by Q/ laser, and hence the probe absorption eventually reaches zero, as seen 
in figure 2. 

iii: Q.I ^ Q^: Dressed states due to combined field: In this condition it is hard to obtain 
analytical expression for the dressed states since the two fields are detuned from each other. 
Only when both fields are on resonance with the atomic transition and hence (0/ — (O^ = 0, we 
can perform rotational transformation which eliminates time dependency in Hamiltonian and 
hence analytically obtain eigenvalues 



Xi =0.0, (19) 

'ki =Ai, (20) 



^3+ = X4+ = , (21) 



A/ 4 


V^' 


+ 2|a/| 


2+^,1 


2 






2 






A/- 


V^' 


+ 2|^/| 


2+^p| 


2 



X3- = ?i4- = . (22) 

We can show that there exists a three-component CPT state for this configuration, given 
by 

identical to one found in reference lfT3l . The probability amplitude for \go) is seen to depend 
upon the ratio of Rabi frequencies of the two lasers Q.p/Q.i. On basis of this, it can be 
interpreted that the occupation probability of |^o) state under these conditions depend upon the 
relative strength of the two lasers. Intuitively extending this to the situation when CO; — cOp 7^ 0, 
one can assume that the superposition state deviates only little from (l23t . In other words, the 
occupation probability of the state |^o) depends only on the relative intensities of the two 
lasers. 

The mixed states : A study of Tr(p^) shows the distribution of atomic population among the 
different CPT states. Figure 5 shows Tr(p^) as a function of co/ — cOp, at different situations. 
The curves labeled 1 to 4 correspond to Q.p — 10^, 1.0, 4.0 and 6.0. It can be noticed that 
Tr(p^) is close to 1.0 throughout when Q.p = 10^, denoting that all the atomic population 
is populating one pure state. As Q.p is increased, the value of Tr(p^) is exactly equal to 1 .0 
only at large values of |(0/ — (Op|, or when (O/ — oOp is exactly zero. For any values in between, 
Tr(p^) deviates from 1.0 and the magnitude of this deviation depends upon value of Q.p. For 
all the data shown in figure 5, (0/ is maintained at atomic resonance and (Dp is varied. The data 
implies that the atomic system is in a pure state when (Dp is far away from (O/, which in other 
words is the situation of a largely detuned Q.p. In such cases, the resonant Q.i establishes a 
stronger coupling and hence traps the atom in its CPT state. On the other hand, when (Dp is 
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exactly equal to (O/, and both beams are on resonance, the atom is in the three component trap 
state (l23t . Tr(p^) being equal to 1.0 when (0/ — (O^ = 0, confirms this fact. For all values in 
between, Tr(p^) < 1.0 indicates that the atom is in a mixed state, which means that the atomic 
population is distributed among many CPT states. 

In figure 6, we plot ImX"^' as a function of (0/ — (O^. For each subset of graphs that are 
shown, Q.I is fixed and each curve corresponds to different Q.p, as described in figure caption. 
Since the pump beam is kept at resonant, variation of co/ — Odp is equivalent to changing the 
probe detuning. It can immediately be seen that for small values of ^/, the probe absorption 
is a characteristic Lorentzian function of probe detuning, as shown in figure 6 (a). When the 
pump beam becomes stronger, as in figures 6(b) and 6(c), the Rabi splitting of levels due 
to dynamic stark shift and onset of Electromagnetically Induced Transparency is noticed by 
presence of side bands and zero absorption for co/ — (Op = 0.0. 

It is now relevant to compare data shown in figure 6 with those in figure 2. The 
corresponding situation of co/ — cOp = 3.0 in figure 2 is equivalent to a point on the abscissa of 
figure 6. Increase of Q.p, as in data of figure 2, is identical to a vertical line in figure 6 over 
which the point of investigation moves. It can immediately be noticed that as ^p is increased, 
the peaks of Im x^^^ curve traverse across this vertical line. Since this peak shape and height 
also change across the curves of figure 6, it is evident that the nonlinear Im j} "> seen in figure 
2 originates in this traverse of peaks across the vertical line. 

In conclusion, we have studied an atomic response to a probe beam which varies 
nonlinearly with respect to the strength of probe beam. The exact trajectory of this variation 
can be controlled using a control beam. The control beam is polarized orthogonally to 
the probe beam and also propagates in a perpendicular direction. This allows a very easy 
control of the two beams without affecting each other. The atomic system, with two sets 
of degenerate triplets is often encountered in many alkali atoms. In short, this is a very 
simple and easy to use system which allows a critical control and nonlinear variation of 
atomic response. Such critical control can have applications in a wide ranging areas of 
information processing, magnetometry etc. A precise control on the nonlinearity of X' "> will 
allow preparation and propagation of soliton pulses within the media as well. Further, the 
phenomena also involves the atom being prepared in different CPT states. This property 
will be very useful in Quantum Computation where a faithful preparation of a state is 
required before performing Computational operations on it. Since these states are CPT 
states and constitute only ground states of the atom, such prepared states are robust against 
intensity fluctuations and decoherence. Thus the state remains little affected till the Quantum 
Computation operations are preformed. Preparations of such states also do not require precise 
pulse shaping techniques. We shall follow up with a work describing such state preparation 
in detail. 
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APPENDIX 

The 36-equations of motion of the density matrix elements are obtained using Liouville 
equation given in the text along with the relaxation operator L p defined by 

[^ P]/m = hm Y, '^yikPkk - 1^ Y« + J^ Jhrn ' Plm, 

where l,m^k = g+,go,g-,e+,eo,e-. We arrange the density matrix equation according to the 
35 X 1 column vector involving the elements of the density matrix p, in the form: 



P = [Pg-eo 


PeQg+ 


Pg-g+ 


Pg-g- Pg+g+ Peog- Peog+ Pg+g- 


Pgoe- Pe+go 


Pe-e+ 


Pe-e- 


Pe+e+ 


Peoeo Pe-go PgQe+ Pe-g- Pe-g+ 


Pe-Ci) Pe+g+ 


Pe+g- 


Peogo 


Pg-gi) 


PgOg+ Pg-e- Peoe- Pe+g+ Pg-e+ 


Peoe+ Pgoeo 


Pgog- 


Pg+gQ. ■ 








Now as explained 


in the text 


, the transformation defined by peg — > pegexp(/(0/f), pg 



Pee and pgg -^ pgg is made. Noting Peg = p*e» we give below the transformed equation 
equivalent to ^ of the text. 



5fPg_eo — [~^' "Yg-eo] Pg-eo -^'^+(Pg-g- " Peo^o) + ^'^P^-Pe-eo " ^'^-Pg-g+ 
5fPeog+ = [^A -y^+eo]Peog+ +'^+Pg-g+ +^'^-(Peoeo -P^+g+) -«^p7+Peoe+ 
5fPg-g+ = ['Ag_g+ -Yg_g+]pg_g+ +/^+peog+ +?^.r_Pe_g+ -/a-pg_eo -«^p7+Pg-e, 
3fPg_g_ = - (2Yg+g_ +2Yg„g_)Pg_g_ +2Yg_e_Pe_e_ +2yg_eoPeoeo +2ygog_Pgogo 

+ 2Yg_g+Pg+g+ +i^+Peog- +i^pT-pe-g--i^+Pg-e,,-i^pT+pg_e^ 

"fPg+g+ = ~ (2Yg-g+ +2Yg()g^)Pg^g^ +2Yg^e+Pe+e+ +2Yg^eoPeoeo +2Yg_g+Pg_g_ 

+ 2Yg()g+ Pgogo +^"-Pe()g+ +'^^pPe+g+ ~'"-Pg+eo ~ '"pPg+e+ 
^fPeog- = V^l -lg-eo]PeQg- + i^+{.Pg-g- - Peoe^■) - i^pPeae- +i^-Pg+g^ 
5fPgoe_ = [/A/ - Ygoe_]pgoe_ " ?a_(pe_e_ - Pg„go) - '^pPgog- + i^-Pe+e^ 
^tPg^,e+ = [/A/ - Ygoe+]pgoe+ - '^-(PgogO - Pe+e+)+ i^+Pe-e+ " '^pPgog+ 
3fPe_e+ = - Ye-e+Pe-e+ + i^pPg-e+ + ?^-pgoe+ - /^pPe_g+ " ?^+pe_go 
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+ iQ.ppe-g--i^-Pe-go 
OtPe+e+ = — {^yg+e+ -'r^ygoe+)Pe+e+ +2Ye_e^Pe_e_ +2Yg^eyPgQeQ +lilpPg^e^ 

+ i^+PgQe+ - i^pPe+g+ - i^+Pe+go 
"tPgogo = ~ {^Jg+go + ^Jg-gQjPgOgO + ■^Ig-goPgOgO + ■^'Yg+goPg+g+ + 2YgQe_pe_e_ 

+ 2Ygoe+Pe+e+ + «^-Pe_go + ^^+Pe+gQ - «^-Pgo<?- ^ i^+Pgoe+ 
dtPe-g- = [iAi-ye-g-]pe-g--i^p{pe-e--pg-g-)+i^-Pgog--i^+Pe-eo 
dtPe^g+ = [i^l - Ye_g+]pe_g+ + i^pPg-g+ + i^-Pgog+ ' i^-Pe^e,, " i^pPe+C- 
dtpe^eo = [«A/ - Ye-eo]Pe_eo + ^^P Pg-^o + ^^-Pgoeo " ''^+Pe-g- - ''^-Pe-g+ 
5fPe+g+ = [«Az - Ye+g+]Pe+g+ -«^p(Pe+s+ - Pg+g+) + /^+Pgog+ -i^-^e+eo 
dtpe+g- = [-i^l - (Ye+ + Yg+)]P^+g+ - i^+Pe+eo + i^-Pgog- + ^^pPg+g- " i^pPe+e^ 
0?Pe+eo = [~'Az — Ye+eoJPe+eo ~ ''^+Pe+g- ~ ''^-Pe+g+ + ''^+Pgoe{) + '^e+g+Pg+eo 
^fPeogO = [-^'Az-Yeogo]Pe()gO + ^'^+Pg-gO+''^-Pg+go ~^'^-Peoe- " ''^+Peoe+ 
^'Pg-gO = ['A - Yg-go]Pg-go + '■^+Peogo + ''^pPe-go - '^-Pg-e_ - i^+Pg-e+ 
^tPg+gQ = [i^l - Yg+go]Pg+go + ''^-Peogo + ''^pPe+go - i^+Pg+e+ - i^-Pg+e- 

and rest of the equations are obtained by, 

Pmn = Pnm 

Where, T± — ± exp{iQOit) . jij is the coherence dephasing term of pjj. Detuning is defined 
as A/ = (Og — cOg — CO/, where ate and oOg are energy of levels \e) and |^) in angular units, pgogg 
is eliminated using equation Q of the text. Then the terms containing r+ form part of the 35 
X 35 matrix M+ .Those containing r_ form part of 35 x 35 matrix M . Terms independent of 
these form the 35 x 35 matrix A. The nonzero elements of /q, however, are given by 



^o(l)='Kog- 


; /o(2) = -/^,„g^ 


/o(4) = /o(5) = Y^og- 


; k{g+) = -i^eog 


W) = i^log. 


; /o(12)=Y.o.- 


/o(13) 


~ Teoe+ 



Then the matrices A, M-^, M and Iq are obtained accordingly. A, M+ and M are 35 
X 35 matrices each and / — is a 35 x 1 matrix in the same order. The nonzero elements of 
A, A/_|_, M are too many in number to be elaborated here. 
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Figure Captions 

Fig. 1 The six level energy level model used in the text. 

Fig. 2 Real and imaginary X'^' as a function of Q.p, for different values of Q/. (a) and (c) 
correspond to absorption and (b) and (d) correspond to dispersion determined from 
equation (14) of the text. The labeled curves are for the Q/ values (1)10^^, (2) 1.0, 
(3) 2.0 (4) 3.0 (5) 4.0. A; = 0.0 and (o, - (o^ = 3.25. y^ = lO^"^) for curves in (a) and (b) 
and Yf = 0. 1 for curves in (c) and (d). All numerical values are normalised to Yeg. 

Fig. 3 Bare states ( — lines), dressed states formed due to Q/ ( — lines) and allowed transitions 
between them (represented by oblique — >). 

Fig. 4 Bare states ( lines), dressed states formed due to strong cOp beam ( — lines) and 

allowed transitions between them brought about by the CO/ beam (oblique — >). 

Fig. 5 Time averaged Tr(p^) for different values of Rabi frequencies as a function of 
relative detuning ((0/ — (Op). The labeled curves correspond to different Q^ values: 
(1)10-^(2)1.0, (3)4.0 and (4)6.0. 

Fig. 6 Absorption of probe beam as a function of relative detuning ((0/ — (Up) for a fixed ^/ 
value and different values of Q^. For figure (a), Q-i = 10^^ and the labelled curves are 
for values of Qp values equal to (1) 10"^ (2) 1.0, (3) 0.2 and (4) 0.6. For figure (b), 
^Z = 1.0 and the labelled curves are for Q.p = (1) 10^^, (2) 1.0 and (3)2.0. For figure (c), 
^i = 4.0 and the labelled curves are for Qp values (1)10^^ (2) 1.0, (3) 4.0 and (4) 6.0. 
The transparency is observed only when Q/ is strong, jc = 10^^ for all the three cases. 
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